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ABSTRACT 

In this paper we describe the results of the first optical kinematic study of the 
planetary nebula BD+30°3639. This system has a central star of the Wolf-Rayet type 
and is believed to be fairly young. Emission line spectra were obtained at high spectral 
and spatial resolution using the Utrecht echelle spectrometer at the William Herschel 
Telescope. These spectra indicate that the main ionized shell of BD+30°3639 appears 
to be evolving in a more complex way than previously thought. 
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1 INTRODUCTION 



The planetary nebula (henceforth PN) BD+30°3639 
(PN G064.7+05.0) is a compact, young PN with a Wolf- 
Rayet type central star. Its surface brightness is among the 
highest of all PNe and the object has be en intensively stud- 
ied at many wavele ngths; fro m X-rays (Arnaud, Borkowski 
fc Harrington 1996 ), opt i cal ( Harrington et al. 1997 ), infra- 
red dLatter et al. 1995; Sicbenmorgen, Zijlstra fc Kriigej 



1994), millimeter (Bachillor ct al. 1991) to radio (Bryce et 



al. 1997a), to name only a few of the many references. 



It is believed to be young because of the low effective 
temperature of the star (~ 40 000 K) and low expansion 
age of the nebula. The central star is special in that it 
shows a hydrogen-deficient, Wolf-Rayet type spectrum (type 
[WC9]). The stellar s pectr um was analysed by Leuenhagen, 
Hamann & Jeffrey (199(:) who derived a surface temper- 



ature {T2/z) of 42,000 K, a luminosity, logj^y L — 4.71, a 
mass loss rate of 1.3 x 10~^ M© yr~^, and a wind veloc- 
ity of 700 km s~^. L and M are distance dependent, and 
these authors assumed a distance of 2.68 kpc to the central 
star (see below). As usual with winds from WR type stars, 
BD-|-30°3639's wind is rather massive. 

Its distance has been determined using the expansion 
proper motions. Two groups used 6 cm VLA images com- 
bined with the optical expansion velocity from the liter- 
ature to determ ine the dista nce. The reported d istances 
are 2.8±fo kpc ([Masson 1989[), 2.680 ± 0.810 kpc (paiiMT 



Terzian fc Bignell 1993 ), and 1.5 ± 0.4 kpc ( Kawamura & 
Masson 19961). These distances are substantially larger than 



1992). Although the proper motion distances should be more 
reliable they are somewhat uncomfortable from the point of 
view of the evolutionary status of BD-)-30°3639. With a dis- 
tance of 2.68 kpc, its luminosity would be around 50,000 Lq, 
highest among PNe. Such a high luminosity would imply a 
massive (> 1.2 M0) central star, which would evolve on 
very short time scales (< 100 years), something which has 
not been observed. Then again none of the published stellar 
evolution models apply directly to WR-type central stars. 

BD+30°3639 is reported to have some more peculiar 
properties, such as a CO expansion velocity substantially 
higher (50 km s~^) than that measu r ed in any other spe cies 
including H2 ( |Bachiller et al. 199lt ^hupe et al. 1998| ). It 
also shows si gns of both C and O-rich dust, the star itself 
being C-rich ( Waters et al. 1998 ). 

Although BD-f30°3639 has been extremely well stud- 
ied, to date no systematic kinematic study of the classical 
optical nebular lines has been done, perhaps partly because 
of the nebula's small size (6" x 4"). However, the kinematics 
are an important element in understanding the PN. They 
can be used to help in deprojecting the object, are also es- 
sential in converting the observed proper motions into dis- 
tances, and can help in the interpretation of the kinemat- 
ics observed at other wavelengths (such as H2 and CO). In 
this paper we report on the first long slit spectroscopy of 
BD+30°3639. 

The lay-out of the paper is as follows: the observations 
and data analysis are described in Sect. 2, in Sect. 3 we 
discuss and interpret the results, and we sum up our con- 
clusions in Sect. 4. 



the previously derived distances using statisti cal methods, 
which lie mostly in the 600 to 700 pc range (Acker et al 
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Table 1. Observational Data 
run 1 



run 2 



run 3 



offset (arcse 



Figure 1. The HST narrow band Ha image of BD+30°3639 

(Harrington et al. 1997) with the EW and NS sht orientations 
indicated. 



2 OBSERVATIONS AND DATA ANALYSIS 

High dispersion, spatially resolved spectra were obtained 
from BD+30°3639 using the William Herschel Telescope 
(WHT) on La Palma combined with the Utrecht echelle 
spectrometer (UES) and the Tekl CCD detector. The ob- 
servations were made on 1994 August 20 and 21. UES was 
used with the 79 grooves mm~^ grating and the maximum 
slit length of 15", long enough to cover the main shell 
of BD+30°3639. However, the wide inter-order separation 
which permits use of a relatively long slit also means that 
full order coverage is not possible with the 1024 x 1024 pixel 
Tekl CCD. Therefore, each slit position was observed at two 
grating positions, covering either end of the 31°' —45"^ orders 
(4740 A - 7415 A) which include most of the bright, kine- 
matically and diagnostically useful optical emission lines. A 
slit width of 252 /xm (= 1.12") was used; this is matched to 
two detector pixels (each 24 /im x 24 /im), giving a velocity 
resolution of 6.12 km s~^. In the spatial direction (with the 
derotator), the scale at the detector is 14.96" mm~^, so that 
1 pixel = 0.36". The 'seeing' remained consistent at ~ 1" 
throughout the observing period. 

The main ionized structure of this PN is remarkably 
rectangular in appearance with its major and minor axes 
aligned approximately along the EW and NS directions re- 
spectively (Fig. nj). Data were obtained from several slit po- 
sitions aligned NS and EW across the nebula. A Thorium- 
Argon wavelength calibration spectrum was obtained at 
each slit and grating position. Flat-field frames were ob- 
tained using a quartz lamp and a wide slit (= 5") exposure 
of the flux standard star SP 1942+261 was also obtained at 
each grating position. 

When these observations were obtained, the A&G sys- 
tem on the UES was not fully operational and it became ap- 
parent after the data reduction process that we could only be 
absolutely certain of the spatial positions of the three runs 
which were clearly centred on the central star (Table |l|) . 

The advantage of obtaining data from several promi- 
nent emission lines in one exposure was largely cancelled 
out by the loss of data due to positional uncertainty and 
the much longer timescale required for the non-standard re- 
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duction process. These factors should be considered if future 
observations of this type are undertaken. 



2.1 Data reduction 

The initial processing was done using standard IRAF rou- 
tines to debias the data and calibration frames and to com- 
bine multiple exposures using standard cosmic ray rejec- 
tion algorithms. No dark frames were obtained since the 
Tek CCD has a very low dark current. The remainder of 
the processing was done using STARLINK routines. The 
data were flat fielded in the usual way before the individ- 
ual position-velocity (pv) arrays were extracted and wave- 
length calibrated. This part of the reduction process was 
non-standard and is described in Appendix A. Although 
some scattering was apparent in the inter-order separations 
around the bright emission lines, this was not removed since 
the 2-d nature of the spectra made it difficult to generate 
an accurate model for subtraction. The scattering is at a 
relatively low level and does not affect the kinematical in- 
terpretation of the spectra significantly. Finally the spectra 
were fiux calibrated and the spectrum of the central star was 
removed from the pv arrays (Appendix B). 

2.2 The main kinematical features 

Examples of the spatially resolved [O ill] and [N ll] emission 
line profiles obtained after the reduction process from the 
NS and EW slits are given in Fig. GI The wavelength axis on 
each image has been converted to radial velocity offset from 
the systemic velocity Usys which was chosen to correspond 
to a heliocentric velocity of —32 kms~^. The spatial scale 
on each image has been converted to an offset in arcseconds 
from the position of the central star. In each case, the stel- 
lar spectrum has been subtracted from the data frame (Ap- 
pendix B). As noted previously, it became clear when the 
data was analysed that the absolute slit positions of several 
of our observations could not be determined. Although we 
obtained unsaturated [N ll]5754 A data at both slit orienta- 
tions, this line is intrinsically faint and so we prefer to show 
the much brighter [N ll]6548 A pv array obtained in run 3 
in Fig. Ba) rather than the [N ll]5754 A pv array obtained in 
run 2. Although the relative brightnesses of these two emis- 
sion lines depends on the electron density of the emitting 
gas, there is no significant difference between the shapes of 
the two pv arrays. Similarly, we obtained [O ill] 5007 A data 
at both E-W and N-W slit positions but again this line is 



© 1998 RAS, MNRAS 000, |-?? 



Kinematics of BD+30P 3639 3 





lU 



L-^ O ^ ^ ^ 

I 

(ossdjd) 1SS||0 




CN O 



(oQSOJd) 1SS||0 



ra 




(0 



CNJ O 




(39S3JD) ^3Sjp 



(DaSOJo) 19SJJ0 



Figure 2. The two-dimensional position-velocity arrays of the (a)[N ll]6548 A (0-5000), (b)[N ll]5754 A (0-100), (c)[0 lll]4959 A (0-15) 
and (d)[0 in]5007 A (0-36) emission obtained from runs 1 and 3. The scale bar corresponds to the ranges indicated, which are in units 
of mjy per pixel. 
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intrinsically weak and the short exposure time for the N-S 
position (run 2) gave a very noisy profile. Hence we prefer 
to show the [O lll]4959 A pv array obtained with a longer 
exposure time (run 3) in Fig. He). There is no problem in 
the direct comparison of these two lines as the [O ill] 5007 A 
and [O lll]4959 A emission lines occur in the ratio of 3:1. 

The main [N ll] shell of emission is more spatially ex- 
tended but less spectrally extended than the main [O ill] 
shell. The velocity ellipses in both high and low ioniziations 
and from both spatial directions are tilted with respect to 
the spectral axis and indeed the [N ll] velocity ellipses ap- 
pear to be almost open ended, with [O ill] emission emerging 
from these gaps. 

The NS and EW slit orientations were chosen because 
these are the projected minor and major axes of the PN 
(see Fig. |l|). For a cylindrically symmetric nebula, these 
projected axes coincide with the direction of the three- 
dimensional (equatorial) symmetry plane and the (rota- 
tional) symmetry axis. However, the spatio-kinematic line 
emission maps shown in Fig. H indicate that this may not 
be the case. The two slit positions show that the blue shifted 
component of the bright [N ll] emission from the NE part 
of the nebula is absent and similarly the red shifted com- 
ponent is absent from the SE part of the nebula. However, 
high speed blue and red shifted components appear in the 
[O III] emission from the NE and SW regions respectively. 
Similar high speed components were also found in the [N ll] 
pv arrays although these components are much fainter than 
the main [N ll] features shown in Fig. bl Since the spatial off- 
sets of the fastest moving components are close to the centre 
of the nebula, this implies a low inclination angle of a pro- 
late nebula. The implications of this are discussed further in 
Sect. 3. 

An [S ll]6716 A/6731 A ratio map was made from the 
data from the NS slit and is shown in Fig. 0, together 
with the [S ll]6731 A position-velocity array. Note that the 
[S ll]6731 A pv array shows the same open ended velocity 
ellipse structure^ as that seen in the [N ll]6548 A pv array 



shown n Fig 



21). The [S ll]6716 A/6731 A intensity ra- 
tio (at a temperature ot 10 000 K), varies trom ~ 1.45 to 
~ 0.45 over an el ectron density range of ~ 10 — lO^cm""^ 
(Osterbrock 1989). The ratio map shows that there is a 
real variation in electron density around the velocity ellipse, 
with the apparently brightest parts having electron densities 
> 2 X lO" cm"^, dropping to 3 - 7 x 10^ cm"^ at the fainter 
ends. The very faint, higher speed, blue-shifted component 
apparently emanating from the northern opening displays 
an electron density < 1000 cm~'^. 



3 DISCUSSION 

3.1 Nebular morphology and deprojection 

The appearance of the line shapes as presented in the pre- 
vious section have some interesting implications for the in- 
terpretation of the morphology of the nebula. We placed 
the slits along the projected minor and major axes of the 
nebula, which in the case of an inclined, cylindrically sym- 
metric nebula would result in a skewed velocity ellipse along 



N a) 

S 



velocity (-;m/s) 






'M 



iloclty (^ni/s) 



Figure 3. (a) Position velocity array of the [S n]6731 A emission, 
scaled between and 200 mjy/pixel. (b) [S ll]6716 A/6731 A 
ratio map scaled between 0.4 and 0.8. The scale bar corresponds 
to these ranges. 



sidering the HST and VLA/MERLIN images one may ques- 
tion whether BD-|-30°3639 is cylindrically symmetric since 
it looks more like a rounded off rectangle, something which 
appears to be inconsistent wit h any kind of proje c tion of a 
cylindrically sy mmetric shell (Bryce et al. 1997a; Harring- 
ton et al. 1997). However, even with this boxy shape one 



would expect the line shapes to contain some information 
about the orientation of the bubble. Combining the NS and 
EW slits shows that the highest redshifts are found in the 
SW and the highest blueshifts in the NE part of the neb- 
ula. The most straightforward interpretation is then that 
the (rotational) symmetry axis of the nebula runs more or 
less diagonally across it roughly from NE to SW, with the 
NE side approaching us. 

The brightness distribution in the line shapes, as well as 
the electron densities derived from the [S ll]6716 A/6731 A 
intensity ratio show that the fastest parts of the PN also have 
the lowest densities, which is consistent with the picture 
of a prolate nebula with faster, low density polar regions. 
The morphology of BD+30°3639 combined with shape of 
the lines shows that we are seeing the nebula almost pole-on. 
Similar line shapes lead to this conclusion in the case of th e 
Ring Nebula NGC 6720 (iBryce, Balick & Meaburn 19941). 



Further evidence comes from a comparison with the model 
line shapes based on hvdr odyii amic simulations presented 
in Chapter 5 of Mellema (1995) and the simple ell ipsoidal 



the major axis and a symmetric velocity ellipse along the 
minor axis. Our results show this not to be the case. Con- 



model based on the distribution of the radio flux ( Masson 



198!: ), which gave a best fit inclination of 18° (but note that 
he assumed the rotational symmetry axis to run EW). We do 
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not think that the low spatial resolution of the data justifies 
more detailed kinematic modelling than this. 

This 'deprojection' is supported by s ome ot her observa- 
tions. As was pointed out in Bryce et al. (1997a) the smooth 
extinction gradient across the nebula (as derived from com- 
paring high resolu tion o ptical and radio emission, see also 
Harrington et al. (1997)), can also be interpreted as being 



due to inclination. Assuming that most of the dust is lo- 
cated in a halo surrounding the entire optical nebula, or in 
an equatorial torus, and that the near side of the nebula 
is the NE, one would expect higher dust extinction in the 
SW, as is observed. The high extinction clump seen in the 
north-eastern part could then be lying along the rotational 
symmetry axis and be the feature in BD-|-30°3639 which lies 
closest to us. 



Graham et al. (199S) reported the detection of an H2 
blob NE of t he ma in optical nebula which was confirmed by 
Shupe et al. (199S). There is some indication of a structure 
at that position in the HST [S 11] /Ha ratio map p resented 
by Harrington et al. (1997). Sahai & Trauger (1998) present 
the HST data in a way which shows the faint small scale 
structures in the halo of BD-|-30°3639. Here one sees ra- 
dial structures in the SW part of the halo and perhaps a 
smoother structure in the NE part. Although none of these 
observations individually could support the claim that we 
are seeing a prolate nebula almost pole-on with the rota- 
tional symmetry axis running NE-SW, combined this seems 
to be the most hkely interpretation. 

Comparing our new optical kine matic data to the H2 
1-0 S{Q) data in Shupe et al. (1998), we do encounter a 
problem. Basically these authors find the opposite kinematic 
behaviour to that exhibited in our data: the H2 emission to 
the east of the nebula is red-shifted by up to 30 km s~^ with 
a faint component reaching 60 km s~^ and that to the west 
is blue-shifted, again by up to 30 km s~^ from the systemic 
velocity. However, our faint, high speed optical components 
are blue-shifted towards the east and red-shifted towards 
the west (the ionized gas at the eastern and western edges 
of the main optical shell is moving almost tangentially to the 
line-of-sight). We checked our slit orientation carefully and 
communicated with the authors of the other paper and the 
effect seems to be genuine. The explanation could be that 
the H2 lies in an expanding fiattened structure around the 
prolate nebula, in which case it would be possible to see both 
positive and negative velocities along the same position an- 
gle, due to the effects of superposition along the line of sight. 
Such a scheme is illust rated in Fig. M. The H2 map presented 
in Shupe et al. (1998) supports this interpretation since it 



shows a ring-like distribution around the ionized nebula. If 
this is the case the expansion velocity of parts of this ring 
would have to be quite high: the projected expansion veloc- 
ity is about 30 km s~^, which with the estimated inclination 
angle leads to a physical expansion velocity which is 2 to 3 
times higher. Unfortunately, the resolution of the H2 data is 
only 50 km s~^, so it is difficult to compare the molecular 
and ionized kinematics one to one. 



3.2 Ionization stratification 

Already in narrow band images it can be seen that there 
is a difference between the low ([ N 11]) and high ([O ml ) 




Figure 4. A sketch illustrating one possible scenario for the kine- 
matical structure of BD-|-30°3639. Moving eastwards from the 
central star, the observer sees firstly the fast-moving, low density, 
optical emission along the approaching polar axis of the PN, then 
the lower ionization, bright optical emission from the main neb- 
ular shell expanding almost tangentially to the line of sight and 
finally the H2 emission from an outer ring of gas with a larger 
component of velocity along the line of sight. 



the high ionization region is smaller than the low ionization 
region. The kinematic data show an even clearer difference. 
The high ionization lines are spatially less extended, reach 
much higher velocities and show a closed shell morphology, 
whereas the low ionization regions are indicative of an open 
ended shell (see Fig. 2). Combined with the deprojection 
suggested in the previous section this means that the prolate 
nebula has fast lower density polar caps which are faint in 
the low ionization lines, but relatively bright in the high 
ionization lines. 



Masson ( |1989| ) suggested that BD-|-30°3639 might be 
like NGC 7027, but seen pole-on. This was based on his 
simple ellipsoidal shell model. Because of the quite severe 
dust extinction in NGC 7027, it is difficult to compare opti- 
cal narrow band images of the two nebulae. There is however 
a nebula which displays exactly the behaviour described in 
the previous paragraph, namely NGC 40. This nebula has 
an ellipsoidal shape with faint polar caps as seen in low ion- 
ization lines, whereas the high ionization lines show a closed 
sheU. Another similarity between BD-f-30°3639 and NGC 40 
is that both have a WR-type central star, t ype [WC8] and 
[WC9] respectively (Leuenhagen et al. 1996). 

When discussing the image and velocity data of 
NGC 40, Meaburn et al. ([l996|) suggested that the sheUs 



ionization regions of the nebula ( Harrington et al. 1997 ) 
© 1998 RAS, MNRAS 000, §-?? 



seen in [N H] and [O III] might have a different origin, the 
first one swept up by an ionization front, and the second by 
the fast wind from the WR-type central star. Nebulae of thi s 
type were foun d in t he numerical models of Mellema ( 1991: ) 
and Mellema (1997). One expects such a configuration to 
occur early in the nebula's evolution, which is consistent 
with BD-|-30°3639. The models show that one clear differ- 
ence between the ionization and wind driven shells is that 
the first does not show much velocity variation as a function 
of latitude, whereas the latter does. In the case of NGC 40 
the [O III] data were not available to test this hypothesis. 
For BD+30°3639 we clearly observe higher velocities in the 
[O III] lines, supporting such a scheme. However, the size 
difference between the nebulae in the fN II] and \0 III] im- 
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ages is rather marginal. As in the case of NGC 40 this means 
that the two sheU scenario remains speculative. 

3.3 Expansion velocities 

BD+30°3639 is one of the radio brightest PNe and its dis- 
tance has been measured by Kawamura & Masson (1996, 
900 ± 200 pc) and Hajian et al. (1993, 2680 ± 810 pc), both 
using the observed proper motions of the expanding shell 
of ionized gas. In both cases the distance was calculated by 
combining the tangential motion derived from observations 
obtained at different epochs with the radial expansion veloc- 
ity derived from optical emission line spectra. Both papers 
obtained an angular expansion rate of ^ 1.3 mas yr~^ for 
the outer edge of the minor axis, and both papers adopted 
an expansion velocity of 22 km s^^. However, Kawamura 
& Masson (1996) then included a correction for a constant 
flux assumption which explains the discrepancy between the 
final distance values. An expansion velocity of 22 km s~^ 
was originally adopted by Masson (1989), in an earlier at- 
tempt to measure the distance to BD-f30°3639, who made 
the point that good, spatially resolved kinematical data was 
not available at that time. 

Our new data provide much more accurate expansion 
velocities. In Fig. 5, spectra from the two central crossec- 
tions of the position-velocity arrays shown in Fig. Wa) and 
(c) are presented. In each case the spectrum was fitted with 
Gaussian profiles and values for the radial expansion veloc- 
ity observed at the centre of the nebula were derived from 
the centres of the two Gaussian components corresponding 
to the two main emission components. This gave values of 
28 ± 1 kms"^ and 35.5 ± 1 kms"^ for the [N ll] and [O ill] 
profiles respectively. The value for the [O ill] expansion is 
considerably higher than the value of 23 km s~^ obtained 
by Sabbadin (1984). The explanation for this discrepancy 
probably lies in the fact that if such a measurement is made 
using lower spatial resolution (i.e. a wider slit or binning 
along the slit), the brighter, but radially slower moving com- 
ponents dominate the spectrum. AUer & Hyung (1995) ob- 
tained a value of ~ 26.8 kms~^ using a 1.2" x 4" slit. In a 
simple proper motion measurement, the calculated distance 
to the object scales directly with expansion velocity, indicat- 
ing that the distance measurements of Kawamura & Masson 
(1996) and Hajian et al. (1993) may have to be revised up- 
wards. This would make the luminosity problem described 
in Section 1 even worse. 

However, a serious problem for these types of distance 
calculation is that they are also dependent on a model which 
relates the radial expansion velocity to the observed tan- 
gential proper motions. If we are observing a prolate nebula 
almost pole-on, the discrepancy between the radial polar 
velocity and tangential equatorial proper motions is at its 
largest. This effectively means that the distances derived be- 
come upper limits. If the outer edge of the PN is actually an 
expanding ionization front (as suggested in Section 3.2), the 
relation between the radial gas velocity and the tangential 
proper motions becomes even more uncertain. 

3.4 Evidence for a high-speed outflow 
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Harrington et al. (1997) presented several emission line ratio 
maps obtained from their WFPC2 images of BD-|-30°3639. 



Figure 5. Spectra obtained by adding together the two central 
crossections from the NS slit for the [O III] 4959 A (upper) and 
[N n]6548 A (lower) emission lines. 



As mentioned above they discovered a compact 'east blob', 
particularly prominent in their [S ll]/Ha ratio map which 
is about twice as far from the centre of the nebula as the 
eastern limb of the main shell and coincides with the region 
of bright H2 emission discovered by Graham et al. (1993). 
Our data contains an interesting feature which could be 
associated with this 'east blob'. The [N ll]6584 A emission 
line from the EW slit is shown in Fig. Ha, on a logarithmic 
scale chosen to show up the very faintest emission surround- 
ing the main, bright, nebular shell (cf. Fig. 0b). The bright 
emission in this line from the main nebular shell saturated 
the CCD chip, but a faint feature extending towards the bot- 
tom left corner can be discerned. This feature lies at about 
5-6" east of the central star and has a velocity ranging from 
~ —60 to —90 kms~^, much larger than the velocities ex- 
hibited by the main nebular shell. Also shown are data from 
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main shells of more evolved PNe (e.g. FLIERS - Balick et 
al. 1993, BRETS - Lopez, Meaburn & Palmer 1993) and 
recently very high spe ed features were obs erved around the 
young PN MyCn 18 (Bryce et al. 19971). In most cases, 
such knots appear in pairs, on opposite sides of the cen- 
tral star, and in cases where more than one pair is observed 
there is often a point symmetry in the pattern. The feature 
in BD+30°3639 is only clearly seen on one side of the cen- 
tral star. Further observations are needed to establish the 
properties of this high velocity feature and to search for a 
possible counter-feature on the opposite side of the nebula. 
In Bachiller et al. (1991;1992) it was noted that the CO 
outflow had a higher velocity than the optical outflow. Our 
data show this to be no longer true. The velocity of the fea- 
ture discussed above reaches 90 km s~^, and even in the 
main nebula these kinds of velocities are reached (see for 
example Fig. 2c). As noted above, the previously reported 
lower optical expansion velocities were based on measure- 
ments of the brighter parts of the nebula. The high velocity 
optical emission does seem to be strongly localised in the NE 
and SW. Bachiller et al. (1992) suggested two explanations 
for the observed CO spectrum: an expanding shell or two 
discrete bipolar blobs. They rejected the latter interpreta- 
tion, but in view of our new results this should perhaps be 
reconsidered. 
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Figure 6. a) The p-v array of the [N ll]6584 A emission from 
the EW slit (run 1), showing a faint, accelerating flow to the east 
of the main nebular shell, b) The p-v array of the [S n]6731 A 
emission obtained from an EW slit position offset slightly north of 
the central star; the accelerating flow is more prominent, c) The 
p-v array of the [N ll] 6548 A emission from the same slit position 
as b). 



another EW slit position, which was slightly offset (probably 
north) from the central star. The fast-moving feature ap- 
pears very clearly in the [S ll]6731 A emission line (Fig. Kp) 
and also, though less prominently, in the [N ll]6548 A emis- 
sion line (Fig. oc) . The intensity observed in the fast moving 
feature expressed as a fraction of the intensity in the bright- 
est parts of the main shell was found to be approximately 
5 times higher in [S ll] than in [N ll]6548 A which is in- 
dicative of shock excitation. The [S ii]6716 A/6731 A in the 
accelerating flow indicates an electron density < 850 cm~^. 
Fast, low ionization features are often found outside the 



4 CONCLUSIONS 

New, spatially resolved, high spectral resolution, position- 
velocity maps have been obtained in both high and low ion- 
ization emission lines from the young PN BD-|-30°3639. The 
low ionization data suggest a prolate, nearly axially symmet- 
ric, open-ended shape for the main nebular shell, with the 
rotational symmetry axis at about PA 30°- 60° in the plane 
of the sky and an inclination of about 20° . The higher ion- 
ization material shows a closed shell structure of a somewhat 
smaller size. This structure may be similar to that seen in 
NGC 40. The high inclination makes it difficult to make the 
connection with the previously published image data, which 
show a remar kably rectangular shap e for the main shell in 
both opti cal (Harrington et al. 1997 ) and radio continuum 
emission (Bryce et al. 1997a). 

The velocity signature for the ionized material se ems to 
be exactly o pposite to that of the molecular material ( shupe 
et al. 1998): the eastern side is mostly blue-shifted in the 



ionized lines and red-shifted in the molecular lines. This be- 
haviour can be explained if the molecular material sits in a 
flattened structure around the prolate nebula. 

New estimates for the expansion velocity determined 
along the line of sight to the central star have been obtained, 
28 ± 1 kms"^ and 35.5 ± 1 kms"^ for the [N ll] and [O ill] 
profiles respectively. These are higher velocities than previ- 
ously supposed, probably because previous observations had 
a lower spatial resolution and therefore incorporated bright 
emission from the shell limbs which has a lower radial veloc- 
ity. These new values can be of use in attempts to obtain the 
distance to BD-|-3Q°3639 using proper motion measurements 
(Kawamura & Masson 1996; Hajian et al. 1993). However, 
it should be kept in mind that the radial expansion velocity 
must be related to the tangential proper motions via a model 
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of the nebular structure and as noted above, our results in- 
dicate that a simple geometry may not be appropriate. 

A fast moving (—90 kms~^), faint feature was detected 
to the east of the main nebular shell. This feature may co- 
incide with structure in th e halo of BD-|-30°3639 detected 
by Harringto n et a l. (1997) in WFPC2 images and by Gra- 
ham et al. (1993) from H2 observations. This feature and 



the high velocity parts of the main shell have velocities close 
t o twice as high as 1: he reported CO velocity of 50 kms~^ 
( iBachiller et al. 199l|) . 
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APPENDIX A: EXTRACTION AND 
CALIBRATION OF INDVIDUAL ECHELLE 
ORDERS 

Neither STARLINK nor IRAF contain specific packages for 
reducing cross-dispersed echelle spectra in the case where 
the 2-d structure of the emission line profiles are to be pre- 
served. This is because, spatially extended, cross-dispersed 
spectra are usually obtained for point-like sources and the 
extended slit thus permits accurate sky subtraction, a pro- 
cess which is folded into the standard reduction packages 
which then produce a 1-d spectrum of the target. The first 
problem with cross-dispersed echelle data is to straighten 
the orders. Each order appears as an approximately verti- 
cal stripe in the raw CCD frame but each 'stripe' is at a 
slight angle to the next one. It is desirable to rebin the data 
such that the 'stripes' become exactly parallel. The largest 
angle between an order and the y-axis was 0.87°, equiva- 
lent to a shift of 15.5 pixels along the x-axis from the first 
crossection to the last. The amount of rotation required to 
straighten the orders was found to be too small for standard 
rotation routines which did not handle the rebinning of data 
sufficiently accurately. The orders were straightened using 
routines which are usually used to SCRUNCH (wavelength 
calibrate) data. We used the continuum spectrum from the 
standard star as a pseudo-calibration arc. This continuum 
spectrum traces out a set of approximately vertical lines. 
Each 'arc line' was identified with the horizontal position 
('rest wavelength') of the continuum spectrum in that order 
in the vertically central crossection of the CCD frame, us- 
ing the ARC2D routine from the TWODSPEC package in 
STARLINK. A dispersion calibration file was then produced. 
This file was used to 'SCRUNCH' the other data files, result- 
ing in a set of CCD frames containing parallel echelle orders. 
The CCDEDIT and TRANNDF routines from the CCD- 
PACK package were used to extract and rotate (through 
90°) each order of interest such that wavelength increases 
from left to right, resulting in a set of pv arrays, each with 
the central star at the same spatial position. 

The PN data frames and standard star frames for each 
order were then wavelength calibrated, following the normal 
wavelength calibration procedure for longslit spectra, using 
the Th-Ar calibration spectra and the routines ARC2D and 
ISCRUNCH. Finally, each PN data frame was corrected for 
air-mass and flux calibrated using standard FIGARO rou- 
tines. The fiux calibration is not entirely rigorous for a vari- 
ety of reasons, primarily because only one wide-slit exposure 
of the flux standard was obtained each night (i.e. one for each 
grating position), although the observing conditions were 
stable. Moreover, the flux calibration is based on an inter- 
polated spectrum derived from calibrated data points spaced 
at intervals of less than one per observed wavelength range 
in a single order. Orders which contain prominent permitted 
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transitions such as Ha and H/3, show wide absorption pro- 
files in the standard star spectrum; these profiles may extend 
further than the limits of the observed wavelength range in 
that order and so may not have been properly corrected for. 
Nevertheless, the internal consistency appears to be good, 
for example the total fluxes in the Ha and H/3 lines, which 
appear in the spectra from both grating positions, were con- 
sistent to within 5% and comparison of the total flux in the 
[N ll]6584 A emission line which appeared at grating posi- 
tion 1 with that in the [N ll]6548 A emission line observed at 
grating position 2 showed that the ratio of fluxes was within 
6% of the theoretical ratio of 3:1. The absolute flux calibra- 
tion was also found to be consistent with other observations. 
An estimate of the total H/3 flux from the nebula, obtained 
by considering the total flux observed through a single slit 
and scaling up from the slit area to the nebular area, is 
8.4 X 10~^* J s~^ m~^. This rather cr ude e stimate compares 



weU with the values in Acker et al. (|l992| ) (9.33 x 10"" J 
s-i m-2) and Harrington et al. ( [1997D (10.5 x 10"" J s"^ 



APPENDIX B: STAR SUBTRACTION 

The central star continuum has been subtracted from the pv 
arrays displayed in this paper, so that the nebular emission 
can be appreciated. This procedure is complicated by the 
fact that the central star of BD+30°3639 is a Wolf-Rayet 
type star with broad emission features superimposed on the 
continuum spectrum. The stellar emission features are about 
an order of magnitude broader than the nebular emission 
line features. The two spatial pixels containing the brightest 
part of the stellar spectrum were flrst extracted from the 
pv array of the order of interest. Any obvious nebular fea- 
tures were clipped out of this spectrum and the resulting 
spectrum was fitted with a polynomial function. The fitted 
function was grown back into a 2-d array. A spatial profile 
of the central star was obtained from the original pv array 
(from a region well away from any nebular emission) and 
this profile was used to modify the regenerated stellar pv 
array in the spatial direction. Finally this pv model for the 
stellar spectrum was normalised to the same fiux scale as 
the original pv array by comparing the flux levels from the 
stellar continuum, again at a position well away from any 
nebular features. This pv array was then subtracted from the 
original array to leave just nebular emission. Careful anal- 
ysis showed that small residuals appeared, corresonding to 
the extended wings of the stellar profiles at positions where 
broad stellar emission features had occured, however, these 
were barely visible over the background noise level and we 
do not believe that they affect our analysis of the nebular 
data. 
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